maximum=max(max(redEr1),max(redEr2),max(redEr3))
minimum=min(min(redEr1),min(redEr2),min(redEr3))

plot, wavenumber1, redEr1, /xlog, /ylog, xr=[1,200], yr=[minimum,maximum], xstyle=1
oplot, wavenumber2, redEr2, color=160
oplot, wavenumber3, redEr3, color='FF0000'x

names=['white -> constant','red -> sqrt(rho)','blue -> rho']
xpos=20
ypos=[8e-04,5e-04,3e-04]

for ii=,2 do xyouts, xpos, ypos[ii], names[ii]
